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Abstract. The origin of highest energy cosmic rays is yet unknown. An appealing possibility 
is the so-called Z-burst scenario, in which a large fraction of these cosmic rays are decay 
products of Z bosons produced in the scattering of ultrahigh energy neutrinos on cosmological 
relic neutrinos. The comparison between the observed and predicted spectra constrains the 
mass of the heaviest neutrino. The required neutrino mass is fairly robust against variations 
of the presently unknown quantities, such as the amount of relic neutrino clustering, the 
universal photon radio background and the extragalactic magnetic field. Considering different 
possibilities for the ordinary cosmic rays the required neutrino masses are determined. In the 
most plausible case that the ordinary cosmic rays are of extragalactic origin and the universal 
radio background is strong enough to suppress high energy photons, the required neutrino mass 
is 0.08 eV < m„ < 0.40 eV . The required ultrahigh energy neutrino flux should be detected 
in the near future by experiments such as AMANDA, RICE or the Pierre Auger Observatory. 



1. Introduction 

The existence of a background gas of free photons and neutrinos is predicted by cosmology. 
The measured cosmic microwave background (CMB) radiation supports the applicability of 
standard cosmology back to photon decoupling which occured approximately one hundred 
thousand years after the big bang. The relic neutrinos have decoupled much earlier, when the 
universe had a temperature of one MeV and the age of just one second. Thus, a measurement 
of the relic neutrinos would provide a new window to the early universe. The predicted 
average number density of relic neutrinos per light (m Vi <C 1 MeV) species i is ~ 56 cm~ 3 . 
This density is comparable with that of photons. However, since neutrinos interact only 
weakly, the relic neutrinos have not yet been detected. 

Recently, an indirect detection possibility for relic neutrinos has been discussed [1-5]. It 
is based on so-called Z-bursts resulting from the resonant annihilation of ultrahigh energy 
cosmic neutrinos (UHECz/s) with relic neutrinos into Z bosons (cf. Fig. 1). The decay 
products of the Z may be detected and yield thereby indirect evidence for the existence of 
relic neutrinos. 

According to recent oscillation measurments (see e.g. [6-9]), neutrinos have non- 
vanishing masses, m Vi . In this case the resonant energy at which the Z production cross- 
section is sizeable reads 
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where M z is the mass of the Z boson. For neutrino masses of 0(1) eV the resonant energy 
is remarkably close to the energies of the highest energy cosmic rays observed at Earth by 
collaborations such as AGASA [10], Fly's Eye [11-13], Haverah Park [14, 15], HiRes [16], 
and Yakutsk [17] (for a review, see Ref. [18]). The possibility that the ultrahigh energy cosmic 
rays (UHECRs) may originate from Z-bursts was first discussed in [4, 5]. 

One of the most outstanding puzzles in cosmic ray physics is the existence of cosmic 
ray events above the Greisen-Zatsepin-Kuzmin (GZK) cutoff [19, 20]. The energy spectrum 
of ultrahigh energy cosmic rays is expected to show a cutoff around 4 • 10 19 eV which is not 
conclusively seen in the data [21]. With the help of the Z-burst scenario one can not only 
detect the relic neutrino background indirectly but also explain the existence of the highest 
energy cosmic rays [22-40]. 

The comparison between the observed and predicted UHECR spectra gives constraints 
on the mass of the heaviest neutrino [32-34, 36]. In this review our quantitative analysis of 
the Z-burst scenario is presented. The requried mass of the heaviest neutrino, as well as the 
necessary UHEz/ flux is determined. 

The organization of this review is as follows. The next section summarizes our present 
knowledge on neutrino masses. Section 3 gives the produced spectrum of protons and photons 
from Z-bursts. The comparison with the observed spectrum is presented in Section 4 and the 
conclusion is given in Section 5. 

2. Neutrino masses 

According to recent experiments, neutrinos almost certainly have non-vanishing masses 
and mixing angles. The evidence for neutrino oscillation is compelling for atmospheric 
neutrinos [6], strong for solar neutrinos [7, 8], and so-far unconfirmed for neutrinos produced 
in the laboratory and studied by e.g. the LSND collaboration [9]. However, neutrino 
oscillations are sensitive only to the mass (squared) differences Am^ = ml. — m\ , not 
to the individual masses, (m U4 >) m V3 > m U2 > m ui , themselves. Only a lower bound on the 
mass of the heaviest neutrino can be derived from these observations, e.g. 
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Figure 1. Illustration of a Z-burst resulting from the resonant annihilation of an ultrahigh 
energy cosmic neutrino on a relic (anti-)neutrino. (Adapted from Ref. [30].) 
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from the atmospheric mass splitting in a three neutrino flavour scenario and 



^Am2 SND > 0.4 eV (3) 

from the still allowed value of the LSND mass splitting in a four flavour scenario. 

For an investigation of the absolute scale of neutrino masses one has to exploit different 
types of experiments, such as the search for mass imprints in the endpoint spectrum of tritium 
beta (J3) decay or the search for neutrinoless double beta (0v2(3) decay. At present, these 
direct kinematical measurements of neutrino masses provide only upper limits, e.g. 



= ^jYlPeil 2 ^ < 2.2-=- 2.5 eV (95%C.L.), (4) 

with U being the leptonic mixing matrix, for the effective mass measured in tritium (3 
decay [41-43], and 



(ITLv 



< 0.33 4 1.35 eV (90%C.L.) (5) 



for the Majorana neutrino mass parameter [44, 45] appearing in Oz/2/3 decay. Recently an 
evidence for neutrinoless double beta decay was reported [46]. The corresponding neutrino 
mass range is: 

(m v ) = (0.11 -=-0.56) eV (95%C.L.) (6) 

Combining the experimental constraints from oscillations and from tritium (3 decay, one infers 
upper bounds on the mass of the heaviest neutrino, 

m V3 < ^ + AmL <2.5eV, (7) 
in a three flavour, and 

m V4 < yjm} + Am 2 LSND < 3.8 eV , (8) 

in a four flavour scenario. 

Further information on the absolute scale of neutrino masses can be obtained through 
cosmological and astrophysical considerations. Neutrinos in the 0.1 -=- 1 eV mass range have 
cosmological implications since they represent a non-negligible part of dark matter. This 
gives the opportunity to put upper limits on neutrino masses from cosmology [47-50]. A 
recent galaxy redshift survey gives an upper bound 

]Tm^< 1.8 eV (9) 

i 

on the sum of the neutrino masses [51]. From the spread of arrival times of neutrinos from 
supernova SN 1987A, coupled with the measured neutrino energies, a time-of-flight limit of 
mp < 23 eV can be derived [53,54], which, however, is not competitive with the direct 
limit (4). According to a recent study, leptogenesis requires that all the neutrinos have masses 
smaller than 0.2 eV [55]. 



3. Z-burst spectra 

The calculation of the Z-burst spectra of protons and photons was performed in three steps. 
First, we determined the probability of Z production as a function of the distance from Earth. 
Secondly, using results from collider experiments we derived the energy distribution of the 
produced protons and photons in the laboratory (lab) system. As a last step the energy losses 
due to the propagation of the protons and photons were taken into account. 
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The prediction of the differential proton flux, i.e. the number of protons arriving at Earth 
with energy E per units of energy (E), of area (A), of time (t) and of solid angle (Q), 

W = dS&j- (10) 

from Z-bursts can be summarized as 

OO fimax 

F p \z(E) = X) / & E P J dr (11) 

i o o 

oo oo 

x J dE Ui F u . (E Ui , r) n D% (r) + J dE Pi F v% (E Dt ,r)n„. (r 
.o o 

x <j ViVi {s) Br(Z -> hadrons) ^H(-) A Pp ( r> Ep , E) , 

with the following important building blocks: the UHECz/ fluxes F Vi (E Vi ,r) at the energies 
E Vi (~ -E^ s ) and at the distance r of Z production to Earth, the number density n Vi (r) of the 
relic neutrinos, the Z production cross section a u . D .(s) at centre-of-mass (cm) energy squared 
s = 2 m^. E^., the branching ratio Br(Z — > hadrons), the energy distribution dN p+n /dE p of 
the produced protons (and neutrons) with energy E p , and the probability P p {r, E p ; E) that a 
proton created at a distance r with energy E p arrives at Earth above the threshold energy E. 

A similar expression as Eq. (11) holds for the differential photon flux from Z-bursts. 
One simply has to insert the energy distribution of photons from Z-bursts and the photon 
propagation function. The latter, P 7 (r, E y ; E) has a different meaning than that of the 
proton's. It gives the expected number of detected photons above the threshold energy E 
if one photon started from a distance of r with energy E 1 . 

The building blocks related to Z-production and decay - o ViVi , the hadronic branching 
ratio, and the momentum distributions dNi/dEi - are very well known. The determination of 
the propagation functions Pi - though CPU intensive - can be done using the precisely known 
microwave backround spectrum and the fairly well known radio and infrared background 
spectra. The first two ingredients, the flux of UHECz/s, F v .(E v .,r), and the radial distribution 
of the relic neutrino number density n Vi (r) are, however, much less accurately known. In the 
following all these ingredients are discussed in detail. 

3.1. Z production and decay 

The s-channel Z-exchange annihilation cross section into any fermion anti-fermion (//) pair 
is given by 

^W^ FMftiM ^' <12) 

where s is the cm energy squared, Y z is the total width of the Z boson, and N eS is the effective 
number of annihilation channels, 

N cS {s) = ^20(s-4m 2 f )^-n f (l-8T 3f q f sin 2 9 W + 8 qj sin 4 9 W ) . (13) 

Here the sum is over all fermions with m/ < with charge qf (in units of the proton 

charge), isospin T 3 / (1/2 for u,c and neutrinos; —1/2 for d,s,b and negatively charged 
leptons), and rif = 1(3) for leptons (quarks (q)). With sin 2 9 W = 0.23147(16) for the sin 2 
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of the effective Weinberg angle and Gp — 1.16639(1) x 1CT 5 GeV -2 for the Fermi coupling 
constant [56], formula (12) gives, at the Z-mass, 

a(uiUi -> Z* -> all qq) \ S=M 2 = 314.9 nb , (14) 

a{y . v . ^ Z *^ all //) | S=M| = 455.6 nb , (15) 

with a branching ratio a"(i/jPj — > Z* — > all qq)/(j{viVi — > — > all //) | S=M |= 0.6912, in 
good agreement with the experimental result [56], 

Br(Z -> hadrons) = (69.89 ± 0.07) % . (16) 

Note that the cross-section is sharply peaked at the resonance cm energy squared 
s = M§. Correspondingly, it acts like a 5-f unction in the integration over the energies E Vi 
in Eqn. (11), and we can assume that the UHECz/ fluxes are constant in the relevant energy 
region. Thus, introducing the energy-averaged annihilation cross section [1,5], 

/A & 
a ann (s) = 2nV2G F = 40.4 nb , (17) 

which is the effective cross section for all neutrinos within 1/2 SE^/E 1 ™ = T z /M z = 2.7 % 
of their peak annihilation energy, we can write 



dE Ui F Vi {E Vi ) a Vi9i {s = 2m Vi E v .) ~ E%° F Vi (E T J*) (a ann ) . (18) 

o 

The next ingredient is the energy distribution of the produced protons and photons in 
Z decay. We combined existing published and some improved unpublished data on the 
momentum distribution, 

V p (x) = ^, x = ^, (19) 

of protons (p) (plus antiprotons (p)) in Z decays [57-61], see Fig. 2 (left). The experimental 
data, ranging down to x w 8 ■ 10 -3 , were combined with the predictions from the modified 
leading logarithmic approximation (MLLA) [62] at low x. The p + p multiplicity is (N p ) = 
Jq dxV p (x) = 1.04 ± 0.04 in the hadronic channel [56]. 

In the cm system of the Z production the angular distribution of the hadrons is determined 
by the spin 1/2 of the primary quarks and thus proportional to 1 + w 2 = 1 + cos 2 9 (here 9 
is the angle between the incoming neutrinos and the outgoing hadrons [63]). The energy 
distribution of the produced protons with energy E p entering the Z-burst spectrum (1 1), 

with y = 2E P /E U , is finally obtained after a Lorentz transformation from the cm system to 
the lab system, 



X 



1 — w 



±y - wJy 2 -(l-w 2 )(2m p /M z y 



y 2 -(l-w 2 )(2m p /M z ) 2 



-wy ± ^y 2 - (1 - w 2 )(2m p /M z y 
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Figure 2. Momentum distributions in hadronic Z decays. Left: Combined data from collider 
experiments on proton (plus antiproton) momentum distribution (solid), normalized to (N p ) = 
1.04, and charged pion momentum distributions (dotted), normalized to (N 7T ±) = 16.99. 
Right: Distribution of protons ("p"; solid), photons ("7"; dotted) and electrons ("e"; dashed) 
in the lab system, in which the target relic neutrino is at rest. 



where m p is the proton mass. The first line comes from the angular distribution, the second 
line is the Jacobian and the third one is the momentum distribution at the inverted momentum. 
The scaled energy distribution Q p , as defined in Eq. (20) and given by Eq. (21), is displayed 
in Fig. 2 (right). 

Neutrons produced in Z decays will decay and end up as UHECR protons. They were 
taken into account according to 

Q P+n (y) = (1 + QM , (22) 

where the neutron (n) (plus antineutron (n)) multiplicity, (N n ) = J 2 dy Q n (y), is ~ 4% 
smaller than the proton's. 

Photons are produced in hadronic Z decays via fragmentation into neutral pions, Z — > 
7r° + X — > 27 + X (cf. Fig. 1). The corresponding scaled energy distribution in the lab 
system, defined analogously to Eq. (20), reads 

Q 7 (y) = f 1 dw Q^o (-^-) , (23) 



-1 1 + w VI + w, 

where the scaled energy distribution Q w o(y), with y = 2 E n o/E„, is given by Eq. (21), with 
m p — > m w o and V p — > V n o, the momentum distribution of pions in hadronic Z decay. For 
the latter distribution, we took the measured one of charged pions V n ± from hadronic Z 
decay [57-61] (cf. Fig. 2 (left)), normalized such that (JV 7 ) = J* dy Q^(y) = 20.97 [56]. 

Electrons (and positrons) from hadronic Z decay are also relevant for the development 
of electromagnetic cascades. They stem from decays of secondary charged pions, Z — > 
ir ± + X — > e ± + X (cf. Fig. 1), and their scaled energy distribution in terms of y = 2 E e /E v 
reads 

Q e ±{y) = f dw f dx 1 = (24) 

x [ 1 2y I v. (■'■) 

\xw + Jx 2 + {2m e /m 7r y 
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Figure 3. Left: The intensity spectrum of the diffuse extragalactic photon background at 
redshift z = (solid). Different estimates of the universal radio background (URB) are also 
indicated: high URB (long-dashed) and moderate URB (short-dashed). Right: Photon energy 
attenuation length Iq at z = corresponding to the photon background shown on the left 
(solid). Variations of Iq arising from different assumptions about the URB are also indicated: 
high URB (long-dashed) and moderate URB (short-dashed). The energy attenuation length for 
electrons due to synchrotron radiation for different magnitudes of the extragalactic magnetic 
fields is also shown (long-dashed-short-dashed). 



where V e is the momentum distribution of the electrons in the rest system of the charged pion. 
The energy distribution Q e ± is also displayed in Fig. 2 (right). 

3.2. Propagation of nucleons and photons 

The cosmic microwave background is known to a high accuracy. It plays the key role in the 
determination of the probability P p (r, E p \ E) that a proton created at a distance r with energy 
E p arrives at Earth above the threshold energy E, suggested in Ref. [64] and determined for 
a wide range of parameters in Ref. [65]. The propagation function P p takes into account the 
fact that protons of extragalactic (EG) origin and energies above « 4 ■ 10 19 eV lose a large 
fraction of their energy due to pion and e + e~ production through scattering on the CMB and 
due to their redshift [19, 20]. In our analysis we went, according to 

dz=-{l + z)H(z)dr/c, (25) 
out to distances R max (cf. (11)) corresponding to redshift z max = 2. We used the expression 

H\z) = Hi [Q M (1 + zf + n A ] (26) 

for the relation of the Hubble expansion rate at redshift z to the present one. Uncertainties 
of the latter, H Q = h 100 km/s/Mpc, with h = (0.71 ± 0.07) xj;^ [56], were included. In 
Eq. (26), Qm an d with fl^ + = 1> are the present matter and vacuum energy densities 
in terms of the critical density. As default values we chose fU/ = 0.3 and = 0.7, as favored 
today. Our results turn out to be pretty insensitive to the precise values of the cosmological 
parameters. 

We determined P p (r, E p ; E) in a range of r < 4000 Mpc, 10 18 eV < E p < 10 26 eV, 
and 10 18 eV < E < 10 26 eV, for several fixed values of the cosmological parameters. 
The simulation was carried out in small (10 kpc) steps in r. For each step, the statistical 
energy losses due to pion/e + e~ production and redshift were taken into account [65]. In 
this connection, the advantage of our formulation of the Z-burst spectrum in terms of 
the probability P p (r, E p ; E) becomes evident. We have to determine the latter only once 
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and for all. Without the use of P p (r, E p ; E), we would have to perform a simulation for 
any variation of the input spectrum, notably for any change in the neutrino mass. Since 
P p (r, E p ; E) is of universal usage, we made the corresponding numerical data for the 
probability distribution (—)dP p (r, E p ; E)/dE available for the public via the World-Wide- 
Web URL http://www.desy.de/~uhecr. 

The determination of the photon propagation function P 7 (r, E 7 ;E), entering the photon 
flux prediction, was done as follows. In distinction to the case of the proton propagation 
function, we used here the continuous energy loss (CEL) approximation which largely 
simplifies the work and reduces the required computer resources. In the CEL approximation, 
the energy (and number) of the detected photons is a unique function of the initial energy and 
distance, and statistical fluctuations are neglected. A full simulation of the photon propagation 
function will be the subject of a later work. 

The processes that were taken into account are pair production on the diffuse 
extragalactic photon background (cf. Fig. 3 (left)), double pair production and inverse 
Compton scattering of the produced pairs. We comment also on synchrotron radiation in 
a possible extragalactic magnetic field (EGMF). For the energy attenuation length of the 
photons due to these processes, we exploited the values quoted in Ref. [66] (see also Ref. [67]) 
and the further ones presented in Fig. 3 (right) which incorporate various assumptions about 
the poorely known universal radio (URB) (from Ref. [68]) and infrared (IRB) backgrounds. 
We shall analyse later the dependence of the neutrino mass and other fit parameters on these 
variations. Note, that, in view of the recent URB estimates in Ref. [68], the ones presented in 
Fig. 3 (left)), which are based on Ref. [69], can be referred to as "minimal" URB. 

The computation of the photon propagation function P 7 (r, Ey, E) was carried out in the 
following way. The energy attenuation of photons in the CEL approximation was calculated 
according to 

where l z (E) = (1 + z)~ 3 lo(E(l + z)) is the energy attenuation length at redshift z. The 
number of photons was assumed to be constant at ultrahigh energies >10 18 eV, due to the 
small inelasticities in this energy range. Below, it was increased in a way to maintain energy 
conservation (except for the redshift contribution): 

dr 

ijE) 

The P 7 (r, P 7 ; E) function was then obtained by integration of these equations. In the 
ultrahigh energy region - which is most relevant for us since we performed our fit to the 
cosmic ray data there - the approximation described above gives the photon flux quite reliable, 
while at lower energies it yields an upper bound. 



d7V 7 = -iV 7 TT7 ^ . (28) 



3.3. UHECv fluxes 

The differential fluxes F u . of UHEz/s are unknown. Present experimental upper limits on these 
fluxes are rather poor. The results of the Fly's Eye [70] and Goldstone lunar ultrahigh energy 
(GLUE) [71] neutrino experiments are summarized of Fig. 4. 

What are the theoretical expectations for diffuse UHECz/ fluxes? More or less guaranteed 
are the so-called cosmogenic neutrinos which are produced when ultrahigh energy cosmic 
protons scatter inelastically off the cosmic microwave background radiation in processes such 
as p7 — > A — > nir + , where the produced pions subsequently decay. These fluxes (for 
recent estimates, see Refs. [3,67,72-74]) represent reasonable lower limits, but turn out to be 



Z-Burst Scenario for the Highest Energy Cosmic Rays 



9 



N 




'I III Mill 



I I I t 



10 



18 



10 



9 



10 20 10 al 10 : 
E„ [eV] 



Figure 4. Upper limits on differential neutrino fluxes in the ultrahigh energy regime. Shown 
are experimental upper limits on F Ve + F Pe from Fly's Eye and on J2 a =e u + ^v a ) from 
the Goldstone lunar ultrahigh energy neutrino experiment GLUE, as well as theoretical upper 
limits on F v + from "visible" ("WB") and "hidden" ("MPR") hadronic astrophysical 
sources. 



insufficient for the Z-burst scenario. Recently, theoretical upper limits on the ultrahigh energy 
cosmic neutrino flux have been given in Refs. [75-77]. Per construction, the upper limit from 
"visible" hadronic astrophysical sources, i. e. from those sources which are transparent to 
ultrahigh energy cosmic protons and neutrons, is of the order of the cosmogenic neutrino flux 
and shown in Fig. 4 ("WB"; cf. Refs. [75,76]). Also shown in this figure ("MPR") is the much 
larger upper limit from "hidden" hadronic astrophysical sources, i. e. from those sources from 
which only photons and neutrinos can escape [77]. A recent detailed simulation showed that 
these bounds can be exceeded if the spectrum of the source is hard enough [74]. Even larger 
fluxes at ultrahigh energies may arise if the hadronic astrophysical sources emit photons only 
in the sub-MeV region - thus evading the "MPR" bound in Fig. 4 - or if the neutrinos are 
produced via the decay of superheavy relic particles. 

In this situation of insufficient knowledge, we took the following approach concerning 
the flux of ultrahigh energy cosmic neutrinos, F Vi {E v ^r). It is assumed to have the form 



where z is the redshift and where a characterizes the cosmological source evolution (see 
also Refs. [3,23,37]). The flux at zero redshift, F Ut (E u . } 0) = F Vi {E n ), is left open. For 
hadronic astrophysical sources it is expected to fall off power-like, F U .(E U .) oc ', 7 > 1, 
at high energies. Due to this fact and because of the strong resonance peaks in the Vii>i 
annihilation cross section (12) at the resonance energies (1), the Z-burst rate will be only 
sensitive to the flux at the resonant energy of the heaviest neutrino. Of course, the latter may 
be nearly degenerate with the other neutrino mass eigenstates, m Ui ~ m v , as it is the case 
for m U3 > 0.1 eV in a three flavour scenario. Correspondingly, our fit to the UHECR data is 
sensitive only to 



F n (E 1/i ,r) = F n (E Vi ,0)(l + z) 



a 



(29) 




(30) 



where the sum extends over the number of mass eigenstates which are quasi-degenerate with 
the heaviest neutrino. Note, finally, that, independently of the production mechanism, neutrino 
oscillations result in a uniform F Vi mixture for the different mass eigenstates i. 
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Figure 5. Left: Mass density fluctuation field 8 along the supergalactic plane as obtained from 
peculiar velocity measurements. Shown are contours in intervals of S — 0.2, surface maps on a 
grid of spacing 500 km s" 1 , corresponding to 5 h^ 1 Mpc, with the height proportional to 5, and 
contrast maps. One recognizes some well-known structures in the nearby volume such as the 
Great Attractor at supergalactic coordinates (SGX — -2000 kms -1 , SGY ~ -500 kms -1 ), 
the Perseus-Pisces complex (SGX — 6000 kms -1 , SGY — -1000 kms -1 ), and the large 
void (SGX ~ 2500 kms~\ SGY ~ kms -1 ) in between. Right: Mass density fluctuation 
field obtained from above data, averaged over all directions, for h = 0.71. The overdensities at 
around 20 and 80 Mpc reflect the Great Attractor and the Perseus-Pisces complex, respectively. 



3.4. Neutrino number density 

The dependence of the relic neutrino number density n Vi on the distance r is the last ingredient 
of the Z-burst spectra. 

The main question is whether there is remarkable clustering of the relic neutrinos 
within the local GZK zone of about 50 Mpc. It is known that the density distribution of 
relic neutrinos as hot dark matter follows the total mass distribution; however, with less 
clustering [40,78,79]. To take this into account, the shape of the n Vi {r) distribution was 
varied, for distances below 100 Mpc, between the standard cosmological homogeneous case 
and that of the total mass distribution obtained from peculiar velocity measurements [80] 
(cf. Fig. 5 (left)). These peculiar measurements suggest relative overdensities of at most a 
factor j v — 2 3, depending on the grid spacing (cf. Fig. 5 (right)). A relative overdensity 
/„ = 10 2 -T- 10 4 in our neighbourhood, as it was assumed in earlier investigations of the Z- 
burst hypothesis [4,5,22-24,35], seems unlikely in view of these data. Our quantitative results 
turned out to be rather insensitive to the variations of the overdensities within the considered 
range. For scales larger than 100 Mpc the relic neutrino density was taken according to the 
big bang cosmology prediction, n u . = 56 • (1 + z) 3 cm -3 . 

4. Determination of the required neutrino mass and the necessary UHECz/ flux 

Our comparison with the observed spectrum included published UHECR data of 
AGASA [10], Fly's Eye [11-13], Haverah Park [14, 15], and HiRes [16], as well as 
unpublished one from the World Wide Web pages of the experiments on 17/03/01. Due to 
normalization difficulties we did not use the Yakutsk [17] results. The latest HiRes data [81] 
are not yet included. We shall take into account the fact that above 4 ■ 10 19 eV less than 50 % 
of the cosmic rays can be photons at the 95 % confidence level (C.L.) [82] 
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Figure 6. The available UHECR data with their error bars and the best fits (long-dashed) from 
ordinary cosmic ray protons originating from our local neighbourhood ("halo background"; 
left) and originating from diffuse extragalactic sources ("EG background"; right), respectively. 
For the latter case, the bump at 4 • 10 19 eV represents protons injected at high energies and 
accumulated just above the GZK cutoff due to their energy losses. The predicted fall-off for 
energies above 4 ■ 10 19 eV can be observed. 



As usual, each logarithmic unit between \og(E/eV) = 18 and log(£7/eV) = 26 is 
divided into ten bins. The predicted number of UHECR events in a bin is given by 

N(i) = £ f El+1 dE [F p{hkd (E) + F p(+l)lz (E)] , (31) 

where 8 ~ 8 • 10 16 m 2 - s • sr is the total exposure of the experiments (estimated from the 
highest energy events and the corresponding fluxes) and where Ei = lfj( 18+ */ 10 ) eV is the 
lower bound of the i th energy bin. The first term in Eq. (31), F p | bkd , corresponds to the diffuse 
background of ordinary cosmic rays from unresolved astrophysical sources. Below the GZK 
cutoff, it should have the usual and experimentally observed power-law form [18]. The second 
term represents the sum of the proton and photon spectra, F p \z + F y \z from Z-bursts. 

The separation of the flux into two terms (one from the power-law background and one 
from the Z-burst) is physically well motivated. The power-law part below the GZK cutoff is 
confirmed experimentally. For extragalactic sources, it should suffer from the GZK effect. 
In the Z-burst scenario, cosmic rays are coming from another independent source (Z-bursts), 
too. What we observe is the sum of the two. As the detailed fits in the next section will show, 
the flux from Z-bursts is much smaller in the low energy region than the flux of the power- 
law background. Correspondingly, the low energy part of the spectrum (between 10 18 5 and 
10 19 3 eV) has very little influence on the Z-burst fit parameters, notably on the neutrino mass. 
We investigated three possibilities for the background term F p | bkd . 

The first possibility is to assume that the ordinary cosmic rays originate from our galactic 
halo (or at least from within the GZK zone of about 50 Mpc). This scenario will be referred 
to as "halo background" in the following. We should note, however, that there are no known 
astrophysical sources in the arrival directions of the highest energy events. This fact makes 
this possibility unlikely. In the halo background case the spectrum of ordinary cosmic rays is 
not distorted by propagation effects, so we can assume the usual power-law behaviour [10, 18] 
(cf. Fig. 6 (left)), 

S F p]hkd (E; A, (3) = ^ ' (Halo bk'd) . (32) 

The second, most plausible possibility is to assume that the ordinary cosmic rays are 
protons from uniformly distributed, extragalactic sources - thus we call it "extragalactic 
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background". In view of the observed, practically uniform distribution of the arrival directions 
of UHECRs, this assumption seems to be phenomenologically more realistic than the halo 
background model. 

The extragalactic background suffers of course from GZK attenuation. Correspondingly, 
we take the above power-law (32), A ■ E~@, as an injection spectrum and determine the final 
spectrum using the proton propagation function P p (r, E p ; E), 

POO fffmax „ 

8 F plhkd (E; A,0)= dE p / dr (1 + z{r)f (33) 



o Jo 



a ( E,y H g5(j^g) (EGbk , d) 



leV VleW v ' dE 
The predicted spectrum of the extragactic background protons shows an accumulation at 
around the GZK scale 4 • 10 19 eV and a sharp drop beyond (see Fig. 6 (right)). 

For comparison with recent work on the Z-burst scenario [37], as a third possibility, we 
have done also fits with no background component, F p | bkd = in Eq. (31), and a larger lower 
end, log(£ min /eV) = 19.4 + 20.0. 

The Z-burst part of the predicted flux (31) has several parameters. As it has already been 
mentioned, this flux is only sensitive to the mass of the heaviest neutrino, this will be denoted 
in the following by m u . The Z-burst spectra can be rewritten as 

£ Fi\ Z (E- B, m v ) =B J dE, J dr (1 + z(r)) 3+a S n (r) (34) 
4m u ( Am u Ei\ dP z {r,E f ,E) 

where a is the cosmological evolution parameter, 5 n (r) is the mass density fluctuation field 
(cf. Fig. 5 (right)), normalized to one, and Qi are the boosted momentum distributions from 
hadronic Z decay, normalized to (N p+n ) = 2.04, fori = p, and to (iV 7 ) = 2 (N n o) + (N w ±) = 
37, for % — 7. After fixing a we have two fit parameters left, the mass m u of the heaviest 
neutrino and the overall normalization B, which may be expressed, using Eqs. (18) and (30), 
in terms of the original quantities entering the Z-burst spectra, as 

I = Br(Z - hadrons) i? max (n Ui ) (<x ann > E™ F?° . (35) 

We can see that the required UHEz/ flux can easily be obtained from the B parameter. 
The background term (when it is nonzero) has two parameters, the (3 exponent and the A 
normalization factor. 

The expectation value for the number of events in a bin is given by Eq. (31). To determine 
the most probable value for m Vj we used the maximum likelihood method and minimized [83] 

the X 2 (P,A,B,m v ), 

log(§)=26.0 

X 2 = £ 2[N(i)-N (i) + N (i)\n(N (i)/N(i))}, (36) 

log(§)=18.5 

where N (i) is the total number of observed events in the i th bin. Since the Z-burst scenario 
results in a quite small flux for lower energies, we took the lower bound just below the "ankle": 
E m - m = 10 18 5 eV. Our results are insensitive to the definition of the upper end (the flux is 
extremely small there) for which we chose log(i? ma x/eV) = 26. The uncertainties of the 
measured energies are about 30% which is one bin. Using a Monte Carlo analysis, we took 
these uncertainties into account and included the corresponding variations in our final error 
estimates. 
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Figure 7. The available UHECR data with their error bars and the best fits from Z-bursts, 
for a strong UHE7 attenuation such that the Z-burst photons can be neglected (a = 0, h = 
0.71, Qm = 0.3,17a = 0.7, z max = 2). Top: The case of an "extragalactic" UHECR 
background. The first bump at 4 • 10 19 eV represents protons produced at high energies and 
accumulated just above the GZK cutoff due to their energy losses. The bump at 2 ■ 10 21 eV is a 
remnant of the Z-burst energy. The long-dashed line shows the contribution of the power- 
law-like spectrum with the GZK effect included. Bottom left : Best fit for the case of a 
halo background (solid line). The bump around 4 • 10 19 eV is mainly due to the Z-burst 
protons (dash-dotted), whereas the almost horizontal contribution (long-dashed) is the first, 
power-law-like term of Eq. (31). Bottom right: The case of no UHECR background above 
log(£ min /eV) = 19.7. 



4.1. Fit results 

Our fitting procedure involves four parameters: ft, A, B and m v . The minimum of the 
X 2 (P,A,B,m u ) function is Xmm a * m ^min which is the most probable value for the mass, 
whereas the 1 a (68%) confidence interval for m v is determined by 

X \f3\ A\ B\ m„) = xl{m v ) = xLn + 1 • (37) 

Here (3' , A', B' are defined in such a way that the x 2 function is minimized in /3, A and B, at 
fixed m v . 

As already mentioned, presently there is no evidence that the observed highest energy 
cosmic rays are photons. Let us start therefore with the assumption (cf. Ref. [32]) that 
the ultrahigh energy photons from Z-bursts can be neglected in the fit in comparison to the 
protons. This is certainly true for a sufficiently large universal radio background, and/or for a 
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Table 1. Results of fits, for a strong UHE7 attenuation (h = 0.71, Am = 0.3, Oa = 
0.7, z max = 2). Top: Assuming a halo UHECR background according to Eq. (32). Middle: 
Assuming an extragalactic UHECR background according to Eq. (34). Bottom: Assuming 
no UHECR background above E min , for different values of the lower end i? m i n of the fit 
(a = 0). 



EG UHECR background + strong UHE7 attenuation 
a m v [eV] xLn A B [3 



-3 


n on +0.20(0.63) 
U - ZU -0.11(0.18) 


25.82 


5.00 • 10 31 


150 


2.465 





n 9n +0.19(0.6l) 
u - zu -0.12(0.18) 


26.41 


5.98 • 10 31 


144 


2.466 


3 


n on +0.19(0.59) 
U - ZU -0.11(0.17) 


26.89 


7.23 • 10 31 


142 


2.467 




Halo UHECR background + strong UHE7 attenuation 




a 


m v [eV] 


V 2 . 
Amin 


A 


B 


P 


-3 


A ,-+2.22(4.80) 
^•^ u -1.64(2.88) 


15.90 


1.46 • 10 43 


1049 


3.110 




3 


q , fi +1.73(4.03) 
J -^ u -1.34(2.32) 
9 .-,+1.45(3.30) 
^• UJ --1.05(1.80) 


15.64 
15.53 


1.62- 10 43 
1.65 • 10 43 


770 

551 


3.111 
3.111 



No UHECR background + strong UHE7 attenuation 



log(£ min /eV) 


m u [eV] 


Xmin 


A 


B 


P 


19.4 
19.5 
19.6 
19.7 
19.8 
19.9 
20.0 


9 98 +0.64(1.46) 
Z - Z °-0.58(1.06) 
, 01 +0.63(1.44) 
^ -0.53(0.80) 
n s .+0.67(1.62) 
U - OO -0.31(0.55) 
n 4n +0.32(0.87) 
u - w -0.16(0.27) 
n ,9+0.41(1.25) 
u,4: -0.18(0.29) 
n 7fi +1.06(2.50) 
u - ' u -0.39(0.58) 
, 7? +1.49(3.47) 
L - 1 1 -1.01(1.47) 


21.81 

16.01 

14.80 

8.03 

7.99 

5.52 

2.68 




1251 
846 
670 
445 
460 
733 
2021 





sufficiently strong extragalactic magnetic field (9(10~ 9 ) G. We shall refer to this scenario in 
the following as "strong" UHE7 attenuation. Our best fits to the observed data for this scenario 
are given by Table 1. We can see that the required neutrino masses are rather insensitive to 
the evolution parameter a. The fits can be seen in Fig. 7, for a = 0. We found a neutrino 

mass of m v = 3.46^J'^2 32) e ^ f° r me " na l° background" scenario, m v = 0.20^' 12(0 is) 
if the background is extragalactic, and m„ = 0.40^[J 'igm 27) ^ there are are no background 
protons above 10 19 7 eV. The first numbers are the 1 a, the numbers in the brackets are the 
2 a errors. This gives an absolute lower bound on the mass of the heaviest neutrino of 
m u > 0.02 eV at the 95% C.L., which is comparable to the one obtained from the atmospheric 
mass splitting in a three flavour scenario, Eq. (2). 

The surprisingly small uncertainties are based on the x 2 analysis described previously. 
The inclusion of the already mentioned 30 % uncertainties in the observed energies by a Monte 
Carlo analysis increases the error bars by about 10%. Note, that the relative errors in the 
extragalactic and in the no background cases are of the same order. This clearly shows that 
the smallness of these errors does not originate from the low energy part of the background 
component. 

The fits are rather good: for 21 non- vanishing bins and 4 fitted parameters they can be 
as low as Xmin = 15.64 and xLm = 25.82 in the halo and the extralactic background case, 
respectively, whereas in the no background case, for E min = 10 19 7 eV, we have 9 bins with 
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Figure 8. The available UHECR data with their error bars and the Xmm ^ 4 fits (solid 
lines) from Z-bursts for an extragalactic UHECR background (long-dashed lines), and a strong 
UHE7 attenuation such that the Z-burst photons can be neglected and only the Z-burst protons 
(dash-dotted lines) have to be taken into account in the fit (a = 0, h = 0.71, Q,m = 0.3, Oa = 



0.7, z n 



2). 



Table 2. Results of fits, for a photon attenuation corresponding to "minimal" URB and 
vanishing EGMF. (a = 0,h= 0.71, Q M = 0.3, fl A = 0.7, z max = 2). 



background 


m v [eV] 


Xmin 


A 


B 


P 


EG 
Halo 


77 +0.48(1.36) 
U - ' ' -0.30(0.51) 
71 +1.40(3.27) 
6 - 1 -1.12(1.96) 


24.52 
15.36 


6.44 ■ 10 31 
1.05 ■ 10 44 


120 
827 


2.467 
3.166 



2 fitted parameters and a xLin = 8.03, see Table 1. In the latter case, however, a remarkable 
dependence of the fitted mass on the value of E min is observed. The ± 2 a fits are shown 
in Fig. 8 for the extragalactic case. According to the expectations the spread is small in the 
region where there are data, whereas it can be quite large in the presently unexplored ultrahigh 
energy regime. 

As the other extreme, let us consider now the case of "minimal" URB on the level of the 
one in [69] and a vanishing EGMF. Then the photons will give a non-negligible contribution 
in the ultrahigh energy region. The corresponding fit results for a = are given by Table 2 
while the best fit for the extragalactic background case is shown by Fig. 9. The values of the 
neutrino mass found in both the halo background scenario, with m u = eV, as 

well as in the extragalactic background scenario, with m„ = 0.77^[j 30(0 51) e V, are compatible 
with the corresponding values found in the case of strong UHE7 attenuation. 

We also performed the fits for other values of the cosmological parameters and relic 
neutrino overdensities. We found that for a wide range (e.g. a = — 3 4- 3, h = 0.61 +• 0.9, 
^max = 2+^5) the results remain well within the errorbars. The only dependence (still within 
the errorbars) is caused by the choice of the background and the type of photon propagation (i. 
e. the strength of the URB and/or EGMF). Fig. 10 contains a summary of the neutrino masses 
for the extragalactic and halo scenarios and three different photon propagations (the already 
mentioned strong UHE7 attenuation and the "minimal" URB and a "moderate" URB case in 
between the two extremes). The overall mass range for the heaviest neutrino in the case of halo 
background is 2.1 eV < m u < 6.7 eV at the 68 % C.L., if we take into account the variations 
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Figure 9. The available UHECR data with their error bars and the best fit from Z-bursts, for 
an energy attenuation corresponding to "minimal" URB and a vanishing EGMF in the case 
of extragalactic background, (a — 0, h = 0.71, fij^ = 0.3,57a = 0.7, z max = 2). The best 
fit (solid line) is the sum of the contributions of the background protons (long-dashed), the 
Z-burst protons (dash-dotted) and the Z-burst photons (short-dashed). 
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Figure 10. Summary of the masses of the heaviest neutrino required in the Z-burst scenario, 
with their 1 a (solid) and 2 a (dotted) error bars, for the case of an extragalactic and a 
halo background of ordinary cosmic rays and for various assumptions about the diffuse 
extragalactic photon background in the radio band (a = 0, h = 0.71, Om = 0.3, ft\ = 
0.7, z max — 2). From left: strong 7 attenuation, moderate and minimal URB. 



between the minimal and moderate URB cases and the strong UHE7 attenuation case. For the 
extragalactic background scenario, the required mass range is 0.08 eV < m u < 1.3 eV. 

Let us consider in more detail the 7 ray spectra from Z-bursts, notably in the ~ 100 GeV 
region. As illustrated in Fig. 11, the EGRET measurements of the diffuse 7 background 
in the energy range between 30 MeV and 100 GeV [84] gives a non-trivial constraints on 
the evolution parameter a. Whereas the high energy spectrum, and thus the neutrino mass, 
is independent of a, at low energies only a<0 seems to be compatible with the EGRET 
measurements, quite independently of different assumptions about the URB. These numerical 
findings are in fairly good agreement with other recent simulations [37]. These photon fluxes 
do not contain any contribution from direct photons emitted by the UHEz/ sources. As they 
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Figure 11. The available UHECR data with their error bars and the best fit from Z-bursts, 
for various cosmological evolution parameters a and an energy attenuation of photons as in 
Fig. 3 (bottom) exploiting a "minimal" URB (h = 0.71, Qm = 0.3, Ha = 0.7, z max = 2). 
Also shown is the diffuse 7 background in the energy range between 30 MeV and 100 GeV as 
measured by EGRET (solid). 



are already close to the EGRET limit, one needs special sources that do not give contribution 
to the EGRET region. 

The necessary UHECz/ flux at the resonant energy E r J s is shown by Fig. 12 together 
with the existing experimental limits and the projected sensitivities of present and future 
observational projets. They appear to be well below present upper limits and are within 
the expected sensitivity of AMANDA [85, 86], RICE [87], and Auger [88]. The fluxes are 
consistent with the ones found in Refs. [37, 38]. 

The astrophysical sources for the required UHE neutrinos should be distributed with 
a <0, accelerate protons up to energies > 10 23 eV, be opaque to primary nucleons and emit 
secondary photons only in the sub-MeV region. It is an interesting question whether such 
challenging conditions can be realized in BL Lac objects, a class of active galactic nuclei for 
which some evidence of zero or negative cosmological evolution has been found (see Ref. [89] 
and references therein) and which were recently discussed as possible sources of the highest 
energy cosmic rays [90-92]. 

It should be stressed that, besides the neutrino mass, the UHECz/ flux at the resonance 
energy is one of the most robust predictions of the Z-burst scenario which can be verified or 
falsified in the near future. 

5. Conclusions 

We have presented a comparison of the predicted spectra from Z-bursts - resulting from the 
resonant annihilation of ultrahigh energy cosmic neutrinos with relic neutrinos - with the 
observed ultrahigh energy cosmic ray spectrum. The mass of the heaviest neutrino should 
be in the range 2.1 eV < m u < 6.7 eV on the 1 a level, if the background of ordinary 
cosmic rays above 10 18 5 eV consists of protons and originates from a region within a distance 
of about 50 Mpc. In the phenomenologically most plausible case that the ordinary cosmic 
rays above 10 18 5 are protons of extragalactic origin the required neutrino mass range is 
0.08 eV < m v < 1.3 eV at the 68 % C.L.. In the case with no background of ordinary 
cosmic rays above, say, 10 19,7 eV we found 0.24 eV < m v < 2.6 eV. The above neutrino 



Z-Burst Scenario for the Highest Energy Cosmic Rays 



18 




I I I I I I L-LLLL 

10 1B 10 19 10 a ° 10 21 10 aa 10 23 



E [eV] 

Figure 12. Neutrino fluxes, F = 3 Yli=i C^U + Fvi)> required by the Z-burst hypothesis 
for the case of a halo and an extragalactic background of ordinary cosmic rays, respectively 
(a = 0, h = 0.71, £Im = 0.3, = 0.7, z max = 2). Shown are the necessary fluxes obtained 
from the fit results of Table 1 for the case of a strong UHE7 attenuation. The horizontal errors 
indicate the 1 a (solid) and 2 a (dotted) uncertainties of the mass determination and the vertical 
errors include also the uncertainty of the Hubble expansion rate. Also shown are upper limits 
from Fly's Eye on F Ve + Fp e and GLUE on J2 a =e uC^* + )> as we U as projected sen- 
sitivities of AMANDA on F V)x + F^ and Auger on F Ua + F 9e . The sensitiviy of RICE is 
comparable to the one of Auger. 



mass ranges include variations in presently unknown quantities, like the amount of neutrino 
clustering, the universal radio background, and the extragalactic magnetic field, within their 
anticipated uncertainties. There is experimental indication that the highest energy events are 
mostly not photons. In this special case, when all ultrahigh energy photons are suppressed 
by a strong enough URB and/or EGMF, the mass range narrows down ( with an extragalactic 
background) to 0.08 eV < m v < 0.40 eV, with a best fit value of m u = 0.20 eV. 

It is remarkable, that the mass ranges required for the Z-burst scenario coincide nearly 
perfectly with the present knowledge about the mass of the heaviest neutrino from oscillations 
and tritium (3 decay , 0.04 eV < m U3 < 2.5 eV, in a three flavour, or 0.4 eV < m U4 < 3.8 eV, 
in a four flavour scenario. 

It is interesting to observe that the recently reported evidence for neutrinoless double 
beta decay [46], with 0.11 eV < (m v ) < 0.56 eV, if true, is compatible with our "favourite" 
extragalactic/strong URB scenario. 

We also determined the necessary ultrahigh energy neutrino flux at the resonance energy. 
It was found to be consistent with present upper limits and detectable in the near future by 
the already operating neutrino telescopes AMANDA and RICE, and by the Pierre Auger 
Observatory. A search at these facilities is the most promising test of the Z-burst hypothesis. 

The required neutrino fluxes are enormous. If such tremendous fluxes of ultrahigh energy 
neutrinos are indeed found, one has to deal with the challenge to explain their origin. It is 
fair to say, that at the moment no convincing astrophysical sources are known which meet 
the requirements for the Z-burst hypothesis, i.e. which have no or a negative cosmological 
evolution, accelerate protons at least up to 10 23 eV, are opaque to primary nucleons and emit 
secondary photons only in the sub-MeV region. 
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